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Abstract. Modern astronomical data on galaxy and cosmological scales have revealed 
powerfully the existence of certain dark sectors of fundamental physics, i.e., existence of particles 
and fields outside the standard models and inaccessible by current experiments. Various 
approaches are taken to modify/extend the standard models. Generic theories introduce 
multiple de-coupled fields A, B, C, each responsible for the effects of DM (cold supersymmetric 
particles), DE (Dark Energy) effect, and MG (Modified Gravity) effect respectively. Some 
theories use adopt vanilla combinations like AB, BC, or CA, and assume A, B, C belong to 
decoupled sectors of physics. MOND-like MG and Cold DM are often taken as antagnising 
frameworks, e.g. in the muddled debate around the Bullet Cluster. Here we argue that these ad 
hoc divisions of sectors miss important clues from the data. The data actually suggest that the 
physics of all dark sectors is likely linked together by a self-interacting oscillating field, which 
governs a chameleon- like dark fluid, appearing as DM, DE and MG in different settings. It is 
timely to consider an interdisciplinary approach across all semantic boundaries of dark sectors, 
treating the dark stress as one identity, hence accounts for several "coincidences" naturally. 



1. Fewer coincidences and dark sectors of the universe 

Dark Matter and Dark Energy are the most fascinating astronomical puzzles presented to 
modern physics of particles and gravitation. To this date, the terms DM and DE are descriptive 
terms without a clear and unique underlying physics. Finding a home for DM and DE in 
the edifice of symmetry-based physics is a challenge. Although analogy is often drawn in the 
literature about DM, DE and Modified Gravity (MG) effects, the more fundamental links of 
these three effects have not been systematically shown. First many MG theories are DE in 
disguise where a special DE field is allowed to be non-uniform and non-minimally coupled 
to the metric. The modification term enters the Einstein's equation on the left or right as 
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known matter, e.g., neutrinos and baryons, G ab is due to the curvature tensor of the metric, and 
M ab is due to extra fields coupled to the metric, and can equivalently viewed as the stress energy 
N ab = ^ G M ab coming from New fields in the DE sectors |2l 131 SI 151 IS1 IZl I&1 IS1 HD1 UH H21 U3| • 
Second the vector, spinor or scalar in MG and DE are special kinds of charge-neutral DM 
fields where DM particles are allowed to self-interact in pairs (e.g., annihilate or be created from 
vacuum), but are forbidden to interact with fields of baryons. This is because for MG and DE 
we have VbM ab = exactly for a Lagrangian construction for the DE or the coupling of MG 
field with the metric. But \7bN ab ~ for DM because of a very small rate of exchange of energy 



when DM and known matter collide. So a DM field with a small cross-section to interact with 
itself and with the baryons is a more general description of MG and DE. 

Third MG, DE, DM fields can be described as a dark fluid with a varying equation of state in 
the framework of General Relativity. The proper treatment involves relativistic hydrodynamics 
rather than CDM-like N-body simulations because DM particles are not point masses moving 
on geodesies, they are packets of energy or fluid which propagates according to the equation of 
motion of the spin field, and the generally non-trivial coupling to itself and to the metric deflects 
the path of DM particles from geodesies. This deflection appears as a fluid-like pressure in the 
stress tensor. Depending on the varying sound speed of the fluid and the Compton wavelength, 
the self-interaction can exhibit as short-range collisions by exchanging massive particles or long- 
range fifth force by exchanging low-energy particle |14[ \TE[ [THJ, [T7\ [T5] or both as in context of 
the lensing and velocity of the Bullet Clusters velocity |19U20[ [2T]. The fluid is a mixture of DE, 
Cold, Warm and Hot DM etc states, and needs not be stable, can make transitions among these 
as the sound speed changes with position and redshift. The energy of the dark fluid is nearly 
conserved, and its cosmic abundance of the dark fluid is nearly fixed at decoupling with the 
radiations. In fact, a plausible candidate for the multi-phased dark fluid is the cosmic mixture 
of multi-flavored neutrinos and antineutrinos, which can make transition among at least several 
energy states, some very low in energy, some very high. Already the uncertain low-energy physics 
of neutrino mass and flavor makes a plausible case for DE in the mass- varying neutrino theory 
(Mass- Varying Neutrinos, MaVaNs, by [22| [23l [2^1 [25] . 

Fourth the self-interaction cross-section must be small and a running function. Many theories 
introducing a zero or a fixed cross-section for self-interaction are excluded by data on either 
small scales (e.g., Cold DM) or large scales (e.g., Collisional DM). Deviations from the ACDM 
model must be small and gradual to preserve its success on Large Scale. Previous multi-field 
models ignore the regularity of DM in galaxies and introduce more degrees of freedoms than 
justified by data. From the perspective of measuring weak lensing, dynamics or SNe distances, 
astronomers detect only the bending of the space-time by the stress energy of all dark fields 
summed together, yet we see a cosmic "coincidence" over 10 30 order in length scale: everywhere 
a characteristic pressure exists T\\ ~ T22 ~ T33 ~ Ac 4 /G ~ (Aml) 2 c 5 /h 3 . This pressure adapts 

to the environment satisfying a Tully-Fisher-Milgrom-McGaugh (TFMM) relation [26l 127} 128]. 
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where the cosmological constant A is related to Milgrom's ao ~ A 1 / 2 ~ tt^—. We expect this 
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small pressure T\\ ~ (100 km/s) 2 x 10 9 M0/kpc 3 can prevent the formation of Cold DM cusps in 
galaxies. A possible home for MG and DE in fundamental physics might well be based on low- 
energy self-interactions of DM, especially what give neutrino mass difference Am 2 ~ (0.03eV) 2 . 

For each galaxy, especially each spiral galaxy, the radial distribution of DM mass Mdm(R) 
and baryons M^ ary {R) could have enjoyed a large scatter as typically seen in CDM simulations 
depending on formation history, but instead the two strictly follow the very tight rule of Tully- 
Fisher-McGaugh globally and the formula of Milgrom |26L [271 [28] at every observable radius 
without significant scatter. These rules hold tightly over many orders of magnitude scatter in 
the gas and stellar surface brightness and a wide range of formation history, including tidal 
dwarfs [33]. 

On bigger scale the universe is in a state of very delicate balance, with an extremely fine- 
tuned ratio of the stress energy coming from the DM and DE, neutrinos and baryons. These 
effects are all within about one order of magnitude of each other at the present instead of some 
100 orders of magnitude apart if they come from truly independent physics. The way to reduce 
fine-tuning is to insist symmetry or restrict degrees of freedom. E.g., these come from a common 
field coupling to all sectors with a cross-section comparable to that of neutrino self-coupling. 

In short we argue for a theory-independent approach to DM, DE and MG. We propose to 
parametrize all these effects and the effects of neutrinos in terms of a Chameleon-like "Dark 
Fluid" {CDF) with an approximately conserved stress energy tensor T a p plus GR. It is a more 
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Figure 1. An ecological view of the Cosmos Most-wanted Pizza (CMP), where different 
dark components (DE, Cold Dark Matter, Neutrino Hot Dark Matter) need not be separately 
conserved. They might be phases of one Dark Fluid. An analogy is drawn with the distribution 
of H2O molecule on earth in phases of water-ice-vapor; any theory with a built-in conservation 
of mass of ice would fail since ice melts. A pure ACDM approach would fail to capture 
the possibility of non-conservation of (Cold) Dark Matter, and while an dynamical/ecological 
approach allows for the possibility of neutrino, CDM and DE to be phase change into each 
other, and explain the coincidence of scales in these three forms of non-baryonic energies in the 
universe. 



systematic approach to DM interactions, since it combines the strength of previous methods 
using purely MOND modified gravity, or Warm DM, or neutrino or Fifth-Force. A natural 
way to extend ACDM is to build in the TFMM relation as the local minimum of the action of 
the self-interacting DM fields. If the DM-DE coupling injects an uneven tiny pressure of order 
Ac 2 /G ~ lO 9 M /kpc 3 x (100km/s) 2 to the random motion pressure of Dark Matter in dwarf 
galaxies, this would suffice to remove the undesired DM cusps in dwarf galaxies. 

The general inspirations of the models are illustrated by Figure. 1, where we seek to address 
the coincidence problems of DE, CDM/MOND and HDM (neutrinos) by one mechanism. In the 
following §2 and §3 we illustrate two possible approaches to unify the treatment of Dark Matter 
and Modified Gravity. §4 illustrates a mathematical equivalence of these two treatments by a 
metric redefinition. 



2. Neutrino approach 

What gives mass to active neutrinos is an open question, and it is not even clear whether 
the mass is constant or not. Interaction with Higgs-like field in other sectors or perhaps self- 
interactions in the neutrino sector (perhaps with sterile neutrinos) give an effective mass m(z, x) 
to the neutrinos observed (chirially left-handed). This argues that the active neutrinos form a 
self-collisional fluid with free-streaming reduced by the pressure. Since the exact mechanism of 
interactions are unknown, one could assume most generally that the effective mass could depend 
on environment. 

The original MaVaNs corresponds to a neutrino fluid of pressure P = —f{n) so that the 
neutrinos can be modeled as a nearly polytropic fluid with an imaginary sound speed in an 
expanding universe. More generally the neutrino pressure P is likely some function of a scalar 
field and its temporal and spatial gradient, P = uEq(3 2 , where (3 = (3(n,n/n, |V$|) depends on 
the density n and its rate of change and spatial gradient (n _1 Vn), and Eq is a constant energy 
scale. 

Unlike MaVaNs, here we assume that certain species of neutrinos have zero mass in the 
absence of perturbations, and the effects of mass appear only in collapsed systems, like galaxies 
and the solar system. We assume that the cosmic fluid of neutrinos has a pressure P due to 
interaction of a current vector nu a of neutrinos and a current vector nu a of anti-neutrinos. 
These two generally do not follow the same equations of motion due to their opposite spin- 
gravity coupling. It is possible for these two currents to interact (creating a pressure) in a way 
reminiscent of how positive and negative charged particles interact in a plasma to keep average 
charge- neutrality on scales larger than the spacing of particles. One expects n ~ n, and a 
coupling like u a u f3 R aj 3 eventually reduces to terms like V a u^V pu a , which in galaxies reduces 
to terms like |V$| 2 . 

To be specific we define (3 to be the sound speed the relativistic neutrino fluid in units of the 
speed of light c, and we adopt a variable sound speed such that 
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where Eq is certain fixed low-energy threshold in the neutrino sector respectively, whose 
microphysics is not specified in detail, but could be thought as certain energy gap in neutrinos; 
note (3 = for the uniform universe with V3> = 0. The quantity 1/(3 carries the meaning of the 
ratio of gravitational wave speed c vs the sound speed c(3 of the neutrino fluid, or the "refraction 
index" of the neutrino fluid, which is assumed to be function of the ratio of the gravitational 
stress 
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vs a characteristic stress nEn = -A^ of the neutrino fluid. One can motivate such a correction 
by saying quantum effects must be considered to model interactions of the metric with a spin 
1/2 neutrino fluid in a box of deBroglie wavelength 1/n ~ l/^o ~ 5mm 3 ; when the gravitational 
energy inside this box is just below or above certain threshold Eq, one might expect that anti- 
particles can disappear or appear (reminiscent of the famous Klein paradox), hence can change 
the energy density of neutrinos and anti-neutrinos, hence the sound speed of the neutrino fluid. 

Interpreting in the picture of the the dielectric analogy of MOND [301 f3Tj one would say 
that the neutrino medium is "polarized" by a variable amount due to a varying gravitational 
stress. The adopted refraction index is such that the perturbations in the stress of the neutrino 
fluid propagate as the speed of light (as in vacuum) in strong gravity, but much slower in weak 
gravity; no propagation or pressure in Minkowski space, where V$ = 0. In the solar system, 
where is much bigger than A, we have (3 = 1. 

We consider a classical Newtonian theory where one minimizes the action S = J dtL given 

by 

S = J dtj dx 3 C, C = p b $ + P g - P, P={nE )(3 2 , (3) 

where the relativistic pressure of neutrino fluid is included by a term nEof3 2 , and P g + p^ 
is the Lagrangian of gravitational field and the baryons of density field pf,, which couples to 
gravitational potential <£. Here we use the non-covariant formulation of gravity as an external 
field, which is the weak perturbation limit of General Relativity. To be clear, R is the scale 
factor, so that in the Lagrangian L we are free to add or drop terms like est J dx 3 R 3 rib and 
est J dx 3 R 3 n, which are the total number of baryons or neutrinos respectively; the est here 
plays the role of the chemical potential. 

The total Lagrangian now resembles the classical MOND of Bekenstein-Milgrom [32J 

C = -^(y-F(y)) + p b $, (4) 

where 

U |V$| 2 

Apply the Lagrangian equation by varying the total Lagrangian with respect to <J>, we obtain 
the MOND-like Poisson equation 

2V [pV<5>\ = 8irGp b , (6) 

where 

(jl = 1 - dF/dy = 1 - (1 + ^)- 3 = 1 - (1 + ^)- 3 . (7) 

It is interesting that a relativistic neutrino fluid with a non-trivial pressure can give the 
physics of MOND. To match with MOND acceleration, we can set 

A = 3a = 3.6 x l(T 10 m/s 2 . (8) 

The resulting p function here matches that of Zhao (2007), which are shown to be compatible 
with solar system data and spiral galaxy rotation curves, especially the Tully-Fisher-Milgrom 
relations. The MOND acceleration scale would corresponds to a mass or energy scale in the 
neutrino sector 
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Eq = —2- ~ 2eV, n = 200cm" 3 . (9) 
8irGn 



Note that any eV range neutrinos with a normal non-relativistic pressure is insignificant to 
galaxy potential. The non-linear coupling of the neutrino pressure to the metric through some 
unspecified quantum effects is the key. In our picture neutrinos are not localized, this can create 
the effects of a ubiquitous dark matter fluid or modified gravity. 

2.1. Relations of neutrinos, Higgs and vector field 

The above example belongs to a more general class of variable mass theory where one minimizes 
the action S = J dtL given by 

S = J dt J dx 3 ^/=g~£, £ = J2 n f m f ~ M p R > ( 10 ) 

where we use the natural units H = 1 and c = 1, and M 2 = is the Planck mass , R is the 
Ricci scalar for the metric g a b, and n/ = *S>*$> is the proper number density of the fermion of 
species / and 



m f = ^g ab f a f b (11) 
is the fermion's effective mass, which is determined by a vector field f a by 

f a = s f (6)m u a + c f {d)u a , ii a = u b V b u a , (12) 

where 9 is an auxiliary field, and u a is a vector field of unit-norm in a metric g ab , defined by 

g ab u a u b = l, g ab = 9abN- 2 , N = — (13) 

m 

i.e., p a = mQU a resembles the four- momentum of a particle of mass thq in the metric g ab , which 
is related to the metric g ab by a conformal transformation factor iV 2 ^> 1. 

To see the relation to the previous neutrino model, we can set N 2 = est in present day 
galaxies, and set sj = est = 0(1), and Cf = for baryons so that baryons have a fixed mass of 
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order Nmo and go on geodesies of g ab , and so that p a = Nmog 00 (1, 0, 0, 0) has a fixed norm 
in galaxies, where goo = 1 + 2$ with $ being the quasi-static gravitational potential. Also set 
Sf = g^(l — 9?) 3 <C 0(1/N) and cj = ip = cos 9 for neutrinos so that neutrinos has a varying 
mass of order m,Q. We can minimize the action against the non-dynamical field 9 to find 

2 

2(c/ -2 + c f ) = (Nm y 2 u a u b g ab ^y= , A = m /N (14) 

in galaxies. This resembles the MOND-like equation in [5j. So cy and the auxiliary field 9 
and the neutrino effective mass rrif = rriQyjF{y) all track y or the acceleration |V<3?| with a 
characteristic scale A = uiq/N. Minimizing against the metric goo or the potential <£, we find 

V- [(/xV$] =AirG^2(n f m f ), fi = l-a^- (15) 

where a = = • These models have properties in between that of classical MOND, 

MaVaNs and Dark Matter. 

In Higgs-doublet models, h and 9 are dynamical fields with a Lagrangian made of following 
kinetic and potential terms 

£h = E( s / /l ) n / + 9 ab V a hV b h + h 2 g ab V a 9V b 9 + V(h) (16) 
/ 



where a potential V(h) ~ cst(h 2 nfmoMp 2 — m^) 2 has a a minimum, can keep a = 1. The 
Higgs field gives mass to baryons and neutrinos through the term Sf(9), but in our case 
the symmetry of the Higgs doublet along the phase angle < < 2n is broken by the term 
Cf(9)ii a spontaneously except in a Minkowski space where u a = 0. The scale h = NrriQ and 
the auxiliary field can take the meaning of the norm and the phase angle of a complex vector 
field Z a = rriQU a exp(i#) and Z a = mou a exp(— iO) with a (dynamical) norm Z a Z h g a \ } = h 2 , as 
discussed in dark fluid models [5]. 



3. The vector approach: a Simple Lagrangian for MOND-like Dark Energy 

Here we illustrate how the roles of both DM and DE could be replaced by a vector field in a 
modified metric theory. This follows from merging two long lines of investigations pursued by 
Kostelecky, Jacobson, Lim and others on consequences of symmetry-breaking in string theory, 
and by Milgrom, Bekenstein, Sanders, Skordis and others driven by astronomical needs [TJ ?.!39j. 

In Einstein's theory of gravity, the slightly bent metrics for a galaxy in an uniform expanding 
background set by the flat FRW cosmology is given by 

9<J> 9^ 

g^dx^dx v = -(1 + -^)d{ct) 2 + (1 - —)a{t) 2 dl 2 (17) 

where dl 2 = {dx 2 + dy 2 + dz 2 ) is the Euclidian distance in cartesian coordinates. In the collapsed 
region of galaxies, the metric is quasi-static with the potential <3?(t, x, y, z) = ^> (t, x, y, z) due to 
DM plus baryon, which all follow the geodesies of g^ u . 

Modified gravity theories are often inspired to preserve the Weak Equivalence Principle, i.e., 
particles or small objects still go on geodesies of above physical metric independent of their 
chemical composition. Unlike in Einstein's theory, the Strong Equivalence Principle and CPT 
can be violated by, e.g., creating a preferred frame using a vector field by, e.g., a unit time- like 
vector field U^ 1 which is designed to couple only to the metric but not matter directly. It doesn't 
violate spatial rotation symmetry since it is time-like. It has a kinetic Lagrangian with linear 
superposition of quadratic co-variant derivatives V(c 2 t/)V(c 2 ?7), where c 2 C/ M is constrained to 
be a time-like four-momentum vector per unit mass by —g^U^U" = 1. The norm condition 
means the vector field introduces up to 3 new degrees of freedom; e.g., a perturbation in the 
FRW metric (Eq{T7|) has c 2 U^ = g^c 2 U" sa (c 2 + <3?, ^p, ^), containing a four- vector made 
of an electric-like potential $ and three new magnetic-like potentials. But for spin-0 mode 
perturbations with a wavenumber vector k, we can approximate — (1,0) ~ (■%, ^-), which 
contains just one degree of freedom, i.e., the flow potential V(t, x, y, z). We expect an initial 
fluctuation of c|k|V ~ |$| ~ c 2 N~ l = 10~ 5 c 2 can be sourced by a standard inflaton; the vector 
field tracks the spectrum of metric perturbation. 

When it comes to writing down a specific Lagrangian of the vector field, simplicity is the 
guide since GR plus simple ACDM largely works. Let's start with forming two pressure terms 



for any four-momentum-like field with a positive norm mc 2 = \J —g a pA a A@ by 
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where the RHSs are co-variant with dimension of acceleration squared, and Vii = A a V a or V Q 
stands for the co- variant derivative with space-time coordinates along the direction of the vector 
A or the dummy index a respectively. From these we can generate two simpler pressure terms 
K and J of the unit vector field U a by 

J = J{U) ~ 0, K = K{U) ~ in galaxies , . 

~ 8nG ' ~ in flat universe 



where the approximations hold for U a with negligible spatial components and nearly flat metric 
(Eq ll7p . Note the J and K are constructed so that we can control time-like Hubble expansion 
and space-like galaxy dynamics separately. The iT-term, with a characteristic pressure scale 

a? 

= Pq in galaxies, is the key for our model. The J-term, meaning critical density, has a 
characteristic scale N 2 Pq ~ 10 10 Po : a t the epoch of recombination z = 1000 when baryons, 
neutrinos, and photons contribute ~ (8, 3, 5) x 10 9 Po respectively to the term J = 3 % G ; so the 
epochs of equality and recombination nearly coincide. 

Now we are ready to construct our total action S = j d 4 x\ — g\^C in physical coordinates, 
where the Lagrangian density 

C = + L m + Lj + L K + {U V U U + l)L m , (20) 

where R is the Ricci scalar, L m is the ordinary matter Lagrangian. For the vector field part, 
L m is the Lagrangian multiplier for the unit norm and we propose the new Lagrangian 

Lj = / J dJAoo(x)| ^ L K = [ K dK\(x)\ ^, (21) 
Jo x 

where the non-negative continuous functions X(x) = (1 + |) 6 — 0, and Aqo = 1 — (l + oo/oo)~ 3 = 
1 — fiB, where pb = 2~ 3 = 1/8. A more fine-tuned parametrization is given in Zhao (2007), 
which passes the BBN constraints better. 

Taking variations of the action with respect to the metric and the vector field, we can derive 
the modified Einstein's equation (EE) and the dynamical equation for the vector field. The 
expressions are generally tedious, but the results simplifies in the perturbation and matter- 
dominated regime that interest us. As anticipated in [30] the ij'-cross-term of EE yields 
\]/ — <J> = for all our models, which means incidentally twice as much deflection for light 
rays as in Newtonian. The tt-equation of Einstein reduces to the simple form 
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, in galaxies (22) 



= fti — j n ma tter-dominated FRW (23) 
3/ifi 3p B 

Here the pressure from the vector field creates new sources for the curvature. The term 
^^llrG^ 1 ^ ™ ^ e Poisson equation acts as if adding DM for quasi-static galaxies. A cosmological 
constant in the Hubble equation is created by 

8^ = - )Kx)d{Pox 2 ) » (3P ) 2 (24) 

For binary stars and the solar system, AirGp— V 2< 3? ~ is true because the gravity at distances 
0.3AU to 30AU from a Sun-like star is much greater than the maximum vector field gradient 
strength ao, so = 0; in fact, |V3>| ~ G ^i G ~ (10 9 — 10 5 )ao, and the typical anomalous 
acceleration is ~ 10 _10 ao, well-below the current detection limit of 10 _4 ao (Soreno & Jezter 
2006). This might explain why most tests of non-GR effects around binary pulsars, black holes 
and in the solar system yield negative results; Pluto at 40 AU and the Pioneer satellites at 100 
AU might show interesting effects. Extrapolating the analysis of [37], we expect GR-like PPN 
parameters and gravitational wave speeds in the inner solar system. 

Near the edges of galaxies, we recover the non-relativistic theory of Bekenstein & Milgrom 
|32j with a function 

H(x) = 1 - \ n {x) ~ p min + x, if x = < 1. (25) 



Note that p(x) — * x hence rotation curves are asymptotically flat except for a negligible 
correction fj, m i n ~ 10~ 15 . In the intermediate regime x = 1 our function with 1 — X n (x) ~ 
(0.55—0.6) for n = 2—5 respectively. Galaxy rotation curves prefer a relatively sharper transition 
than p{x) = x/(l +x) = 0.5 at x = 1 [31] where we can identify 9b I {.9dm + 9b) = fJ>(x). So our 
model should fit observed rotation curves. 

For the Hubble expansion: the vector field creates cosmological constant-like term ~ 9Po 
below the zero-point of the energy density in the solar system because the zero point of our 
Lagrangian (Eq l2ip is chosen at N 2 Pq < K < +oo. During matter domination, the contribution 
of matter 8nGp and Ao to the Hubble expansion H 2 (Eq l23|) is further scaled-up because the 
effective Gravitational Constant G e jf = G/ps = 8G > G. Coming back to the original 
issue of the 3 : 1 ratio of matter density to our cosmological constant, Eq. (|23p predicts that 

9P^ . 4(1+ z) Pq ^^v, ; + „ + v, Q ^l^o^^ q ■ n _l ^3 
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which is close to the desired 3 : (1 + z) 3 ratio. Adding neutrinos 



makes the explanation slightly poorer. So the DE scale is traced back to a separate coincidence 
of scale, i.e., the present baryon energy density pbC 2 ~ 4Pq, where Po contains a scale ao for the 
anomalous accelerations on galactic scale. This model predicts that DE is due to a constant of 
vacuum. 

In our model, the effective DM (the dog) follows the baryons (the tail) throughout the 
universal (1 + z) 3 expansion with a ratio set by the parameter \ib- To fit the ACDM-like 
expansion exactly, we note the Hubble equation for a flat FRW cosmology with vector field and 

tt h h 2 ^ Q v h 2 0.07eV ^ Q ph h 2 
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The 2nd equality fixes /z^ 1 = (8 — 8.4) if we adopt (iq/c Hq/6 ps 12km/s/Mpc and 
= (0-25 — 0.3). The 1st equality would predict an uncertain but very small neutrino 
mass m u ~ ±0.3eV, consistent with zero. So the role of neutrinos in uniform expansion can be 
completely replaced by the vector field. 



4. A unified framework for Dark Matter, Dark Energy and Modified Gravity 

Finally we illustrate the relation between Modified Gravity and Interacting Dark Matter. 



4-1- Einsteinian gravity with an interacting Dark Matter field 

Let's consider Einsteinian gravity but with the normal matter (of standard model of particle 
physics) being coupled to the field of a dark matter particle. Let the dark particle be spin-1, 
hence it is described by a vector field, hence with 4 degrees of freedom. Unlike the spin-1 photon 
field, which is a massless gauge-invariant vector field with zero expectation value in vacuum, the 
dark particle vector field is given a unit norm, hence it is a massive field and has a non-zero 
expectation value in vacuum. The vector field is given self-coupling and coupling to normal 
matter, which break the gauge-invariance. 

The gravity sector is now simply described by a metric g a b with a sign convention (+,-,-,-) 
and its associated Ricci scalar R, hence the Einstein-Hilbert action plus matter, 



S — Sn ~\~ Sn 



dx q 
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R 
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+ L r 



Lm — La + Li n t + Lj 



(27) 



where matter is consisted of the Lagrangian density Lj for a pressure-less matter fluid with 
matter flux vector J a = pu a , and La for the dark matter vector field A a , and an interaction 
term of the two vector fields. Specifically 



Lj = ^g ab J a J b -0V a J a , 



(28) 



where <fr is a Lagrangian multiplier field for the conservation of matter flux J a = pu a of a 
collisionless dust of density p and four- velocity u a . It interacts with the vector field A a via 



L int = C 2 VBC^/g ab J^ gcd A^Ad - (1 - B)( to JM^ - yJ 9ab J a J b (29) 
where B and C are coupling constants. The vector field contributes via 



m 



F ab F ab + (l-g ab A a A b )X 



ab 



g bc V a A c - g ac V b A c , 



(30) 



which consists of a kinetic photon- like term from Faraday tensor F ab , and a massive potential 
in the form of a unit norm constraint for A a and a Lagrangian multiplier A. And to make the 
argument simpler we assume the mass m and A oc m 4 is very small, f3 = lQirGm 2 <C 1, so that 
we can neglect the dark matter term La, keeping only the interaction term Lj nt . 

Note that our Lagrangian Lj depends on g ab in a non-linear fashion. However, non-linearity 
is not a sufficient condition for modified gravity. It can be shown that the stress tensor associated 
with Lj is T ab = ~?j=^5p = JaJb{g ah JaJb) 1 / 2 = pu a u b , as expected for a collisionless dust. 

Note that we cannot observe the dark matter A a field directly, it is observed through its 
interaction with baryonic dust J a . Here we considered only a species of dust. If we generalize 
for the interaction/coupling to be the same for all species of the baryonic dust, then we cannot 
detect the dark matter through differential measurements of the baryons (the strong equivalence 
principle). However, the dark matter needs not track the baryonic dust exactly, so we don't 
expect the baryonic mass center to coincide with its kinematic (gravitational) center (Kesden & 
Kamionkowski 2006). 



4-2. "Modified" gravity in redefined metric 
Alternatively one can redefine the metric 



9ab = (gab - (1 - B)A a A b ) C 
g ab = [g ab - (1 - B- 1 )A a A b ) /C. 



(31) 
(32) 



The nice feature is that the matter action is now simplified to that of a pure matter field J a , 
with a new Lagrangian density 



Lj 



^J~9abJ a J b ~ 0V a J a 



(33) 



which is completely decoupled from the vector field A a . However, the new Ricci scalar, formed 
out of second derivatives of the new metric, differs from the old Ricci scalar by a K-term, so 
that the gravity is "modified" with an effective gravitational constant G given by 



G 



'-g 



G C^g 



BC. 



(34) 



Rewriting the action S in the new metric, we find an action resembling that of a "modified" 
gravity, specifically 



where 



S = S g + Sj = - dx A ^g 



K = 



R + K 
T&rrG 



-I 



dx \/-gLj, 



K±V a A' m V b A' n ] + (g ab A' a A' b - 1)A, 



(35) 
(36) 



A' a = A a /VBC (37) 
K% = {d x g ah + c' 4 A' a A' b )g mn + (c' 2 5 a m 5 b n + c' 3 5 a n 5 b m ) (38) 
where A is a Lagrangian multiplier field, 

1 - B- 1 , 

c[ = -c 4 = c' 2 + B(5, -c' 3 - c[ = c ' 2 = B -1. (39) 

where we have considered more generally when the mass m of the vector field is not small, the 
c-parameters are functions of (3 = \QirGm 2 as well [36j . 

The redefined action is that of a special case of the Einstein- Aether modified gravity, where 
there is no interaction between normal matter and the vector field (called aether). Assuming A a , 
J a and g ab as independent freedoms as conventionally done, the Einstein equations are obtained; 
taking its trace we have 

2R - ~g ab T± = 8ttGT j , f J = % = b (40) 

VK^T b 

where T J is the trace of the stress tensor of the collisionless dust, and R is the Ricci scalar, 
— is the part of Einstein tensor involving a fairly lengthy expression of second derivatives of 
the vector field A a and the metric g. For the uniform expansion of the universe in co- moving 
coordinates (t,x,y,z), the metric is given by 

g ab dx a dx b = dt 2 - a{tf{dx 2 + dy 2 + dz 2 ), (41) 

the matter current J a = (p, 0,0,0) and A' a = (1,0,0,0) with a preferred time-like direction. 
The Hubble equation is unchanged, 

3«f = 8^, ±- , g/g - ° (42) 

\a{t)dt) nc Ac 1 + (ci + 3c 2 + c 3 )/2 yfB 

except for a correction [ic of the effective gravitational constant. For weakly perturbed metric 
near a static galaxy or the solar system with 

g ab dx a dx b = (1 + 2<S>)dt 2 - (1 - 2$>){dx 2 + dy 2 + dz 2 ), (43) 

the matter current J a = (p, 0,0,0), and the vector field A' a = (1/y/goo, 0, 0, 0) with a preferred 
time-like direction, and the Poisson equation is changed to 

V 2 (2^$) = SttGp, — = = CVB (44) 

UN l-(ci + c 4 )/2 

with a factor /ijy correction of the effective gravitational constant. 

One can create appropriate amount of dark matter-like effects by selecting the factors pn 
and pc in solar system, in galaxies and in matter-dominated cosmology. The above arguments 
can be generalized to the case where pjy and pc are scalar fields, which varies with redshift and 
position to resemble the V — A model and F{K) models of co-variant MOND [3] and various 
modified gravities. E.g., the modified source gravity of Carroll et al. can be recovered by models 
with pc = pb = C being non-dynamical scalar fields tracking the Ricci scalar through an added 
potential term V(C) in the interaction Lagrangian. 

In short, the above example suggests that dark matter and modified gravity can be the two 
view points of the same phenomena. Our vector field can be viewed as spin-1 dark matter 



field in Einsteinian gravity and g a b, and this dark matter field A a is co-variantly coupled to the 
luminous matter current field J a . However, this vector field is not exactly the cold dark matter, 
and in fact it does not condense in galaxies in our model, but it does decelerates the expansion 
of the universe, and contributes to the critical density. Alternatively, one can view our vector 
field as a field modifying the gravity sector of metric g. In this case, the luminous matter is 
decoupled from the vector field A a . 

In short the coincidences of scales of DM, DE and Neutrinos are intriguing. We advocate 
that it is theoretically satisfying to find a unified solution to these problems at a fundamental 
level. 
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